Electron-Phonon Interactions in the W = Pairing Scenario 
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We investigate the interplay of phonons and correlations in superconducting pairing by introducing 
a model Hamiltonian with on-site repulsion and couplings to several vibration branches having 
the Cu-0 plane of the cuprates as a paradigm. We express the electron-phonon coupling (EP) 
through two force constants for O-Cu and O-O bond stretchings. Without phonons, this reduces 
to the Hubbard Model, and allows purely electronic W — pairing. A W = pair is a two-body 
singlet eigenstate of the Hubbard Hamiltonian, with no double occupancy, which gets bound from 
interactions with background particles. Indeed, this mechanism produces a Kohn-Luttinger-like 
pairing from the Hubbard repulsion, provided that its symmetry is not severely distorted. From 
the many-body theory, a canonical transformation extracts the effective two-body problem, which 
lends itself to numerical analysis in case studies. As a test, we use as a prototype system the 
Cu04 cluster. We show analytically that at weak EP coupling the additive contributions of the 
half-breathing modes reinforce the electronic pairing. At intermediate and strong EP coupling and 
U ~ t, the model behaves in a complex and intriguing way. 
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I. INTRODUCTION 

While the Frohlich mechanism of conventional supercon- 
ductivity is driven by phonon exchange, the pairing mecha- 
nism in highly correlated, narrow-band systems could have 
a predominantly electronic origin [1] and the Cu-0 plane of 
cuprates is the most discussed example. Although this re- 
mains a very controversial issue, most authors probably ac- 
cept at least the conceptual importance of a lattice coun- 
terpart of the Kohn-Luttinger idea [2] that attractive interac- 
tions result from mere repulsion. The Renormalization Group 
approach [3] to the Hubbard Model shows that such a super- 
conducting instability in the d^2_y2 channel is dominant near 
half filling, confirming the results obtained with the FLEX 
approximation [4] . One definition of pairing is A < 0, where 



A(iV + 2) = E{N + 2) + E{N) - 2E{N + 1), 



(1) 



and E{N) is the ground state energy of the system with N 
fermions; this criterion is suitable for finite cluster calcula- 
tions by exact diagonalization methods. A < was indeed 
observed in particular Hubbard clusters [5] [6], but in many 
other examples with on-site repulsion on every site a large 
A > was found [7] [8]. The W = theory [9] [10] [11] gives 
a systematic method for producing and analizing examples of 
singlet pairing by repulsive interactions; it also allows validat- 
ing |A| as the pairing binding energy. In this framework the 
non abelian symmetry group of the underlying graph and the 
resulting degeneracy is crucial for the pairing mechanism. The 
fillings and the symmetry channels where the W = pairing 
can occur are determined in full generality by the W = the- 
orem [12]; these symmetries achieve the same result as high 
angular momentum and parallel spins in the Kohn-Luttinger 
[2] continuum approach. 



Anyhow, a purely electronic theory misses practically and 
conceptually important features of this complicated problem. 
First, many high- To compounds exhibit a quite noticeable 
doping-dependent isotope effect [13], suggesting that electron- 
phonon (EP) interactions are important and should be in- 
cluded in the theory. In addition, there is experimental evi- 
dence [14] that the half breathing Cu-O bond stretching mode 
a.t k — (tt, 0), (0, tt) is significantly coupled with the doped 
holes in the superconducting regime and its contribution may 
be relevant for the d^2_y2 pairing [15] [16] [17]. A radical, yet 
serious criticism of all electronic mechanisms was put forth 
by Mazumdar and coworkers [18]. They suggested that any 
pairing in Hubbard clusters is of doubtful physical interpre- 
tation due to the neglect of the lattice degrees of freedom 
and the Jahn- Teller (JT) effect. They argued that JT distor- 
sions might well cause a larger energy gain of the system with 
A^ + 1 particles, and that could reverse the sign of A obtained 
at fixed nuclei; in this case the pairing would be just an arti- 
fact of the Hubbard model. This even led the authors to the 
conjecture that any A < due to an electronic mechanism is 
just a finite size effect, which vanishes for large systems like 
the JT effect does. Below, we shall show that the Mazumdar 
et al. argument [18], based on a static Jahn- Teller approxima- 
tion, can break down in a peculiar and nontrivial way when 
more flexible wave functions are allowed; thus, phonon pair- 
ing and W = pairing can be compatible, depending on the 
symmetries of both pair and vibration and on frequencies. 

The Hubbard-Holstein model, where electrons are coupled 
to a local Einstein phonon, is a simple way to include both 
strong electronic correlations and EP interactions. Much is 
known about the possibility of a superconducting phase in 
this model. Pao and Schuttler [19] applied the numerical 
renormalization group techniques within the FLEX approxi- 
mation and found that in the square geometry s-wave pairing 



is enhanced by phonons, while d^2_y2 pairing is suppressed. 
In tlie strong EP regime a Lang-Firsov [20] trasformation 
maps the Hubbard-Holstein model in an effective Hamilto- 
nian for hopping polarons with a screened on-site interaction 
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U = t^ ~ 771 where U is the Hubbard repulsion, uo is the 
phonon frequency and g is the EP coupling constant. If over- 
screening is attained, f/ < becomes an effective attraction 
and one gets bipolaronic bound states [21] [22]; however it is 
still unclear whether they can exist as itinerant band states 
[23] [24]. Petrov and Egami [25] found A < in a doped 8- 
site Hubbard-Holstein ring at strong enough EP couling, while 
otherwise the normal repulsion prevails. This situation is un- 
avoidable in 1 dimensional repulsive Hubbard model, where 
no superconducting pairing exist. 

In this paper we take the view that one of the sound ex- 
perimental facts about the CuO plane in all cuprates is that 
the geometry permits W = G pairs which avoid completely 
the strong hole-hole repulsion. It is therefore highly plausi- 
ble that such pairs are important ingredients of the theory, 
provided that there is a way out of the Mazumdar et al. [18] 
argument. Anyhow, it is not obvious that the phonons will 
reinforce the attraction while preserving the symmetry. More 
generally, some vibrations could be pairing and others pair- 
breaking. When lattice effects are introduced in the W = Q 
scenario, the situation is very different from the Petrov and 
Egami [25] model, when, as in the conventional (Frohlich) 
mechanism, phonons overscreen the electron repulsion; what 
happens if electronic screening already leads to pairing? 

To address these problems we use an extension of the Hub- 
bard model in which bond stretchings dictate the couplings to 
the normal modes of the C4„-symmetric configuration. This 
is physically more detailed than the Hubbard-Holstein model, 
and does not restrict to on-site EP coupligs that would be im- 
paired by a strong Hubbard repulsion. 

The plan of the paper is the following. After introducing 
the model Hamiltonian in the next Section, we devote Sect. Ill 
to a detailed derivation of the effective interactions between 
holes in the W = G pair, which is obtained by extending a 
previous Hubbard Model treatment. Our canonical transfor- 
mation approach is quite general for weak EP coupling and 
corresponds to the inclusion of all diagrams involving one- 
phonon and electron-hole pair exchange due to correlations. 
We specialize in Sect. IV to the prototype CUO4 cluster, de- 
scribing electronic states and vibration modes. The effective 
interaction is calculated explicitly in Sect.V. Next, we de- 
velop a theory based on the Jahn- Teller operator in Sect. VI; 
in this way we want to test the reliability of that approxima- 
tion in modeling the behaviour of VF = pairs in the presence 
of Jahn- Teller active modes. The numerical results of the full 
theory are then exposed and discussed in Sect. VII; the exact 
data for realistic vibration frequencies disagree from those of 
Sect. VI but are in accord with the canonical transformation 
approach of Sect. V. The agreement is excellent at weak cou- 
pling, but the analytical approach is qualitatively validated 
also at intermediate coupling. Finally Sect. VIII is devoted to 
the conclusions. 



II. MODEL 

We start from the Hubbard model with on-site interac- 
tion U and expand the hopping integrals tij{ri, Vj) in powers 
of the displacements pi around a C4„-symmetric equilibrium 
configuration 



a 



OTij [Yi , Vj 



drf 



P, + 



(2) 



where a — x,y. Below, we write down the p" in terms of 
the normal modes g,, „: p" ~'^ S"^{i) qr^v, where 77 is the 
label of an irreducible representation (irrep) of the symmetry 
group of the undistorted system and ;/ is a phonon branch. 

Thus, treating the Cu atoms as fixed, for simplicity, one 
can justify an electron-lattice Hamiltonian: 

Hel-latt — Ho + Vtot • (3) 

Here Ho = Hq + Hq is given by 

Ho = 22 ^'^v^i^bl^^br,^,, + 2J tlj{ri,rj){cl^c.ja- + h.c), (4) 

where ujr).ii are the frequencies of the normal modes with cre- 
ation operator foj,_^, while c\^ creates a fermion of spin a in site 

i. Moreover, let M denote the O mass, £„ u = \v\^TTr — 

with A^ V numbers of order unity that modulate the EP cou- 
pling strength. Then, Kot = V + W reads 



Kot = 2j?').''(^ii7,-' + br,,v)H,^,„ + Uy^nj^riiy, 
the Hr^^u Operators are given by 



(5) 



Hr, 



+S^M) 



EE 






(») 



dtiAYi,Yj 



drf 



dUj {r„rj) 



9rf 






+ h.c: 



(6) 



In previous work we have shown that the pure Hubbard 
model Hh = Hq + W defined on the Cu-0 plane and on 
the simple square as well admits two-body singlet eigenstates 
with no double occupancy on lattice sites. We refered to them 
as W = pairs. W = Q pairs are therefore eigenstates of the 
kinetic energy operator Hq and of the Hubbard repulsion W 
with vanishing eigenvalue of the latter. The particles forming 
a Vy = pair have no direct interaction and are the main 
candidates to achieve bound states in purely repulsive Hub- 
bard models [26] [27] [28]. In order to study if the W = 
can actually form bound states in the many-body interact- 
ing problem, we developed a canonical transformation of the 
Hubbard Hamiltonian [26], which enables us to extract the 
effective interaction between the particles forming the pairs. 
Pairing was found in small symmetric clusters and large sys- 
tems as well. 

In the next Section we wish to derive an effective inter- 
action between the particles in the pair suitable for H^i-iatt, 
by generalizing the canonical transformation approach of Ref. 
[26]. 



III. CANONICAL TRANSFORMATION 

In this Section we assume that the system has periodic 
boundary conditions with particle number A'^; we denote the 
phonon vacuum by |0)) and the non-interacting Fermi sphere 
by |$o(A'^))- The creation operator oi a W = pair is ob- 
tained [12] by applying an appropriate projection operator to 
cLcJ., |, where the labels denote Bloch states. If we add a 
W = pair to |$o(Af)) (g) |0)), the two extra particles, by def- 
inition, cannot interact directly (in first-order). Hence their 
effective interaction comes out from virtual electron-hole (e- 
h) excitation and/or phonon exchange and in principle can be 
attractive. To expand the interacting (Ai'+2)-ferniions ground 
state |^o(-^ + 2)), we build a complete set of configurations 
in the subspace with vanishing z spin component, considering 
the vacuum state |$o(A'')) ® |0}) and the set of excitations 
over it. 

We start by creating W — pairs of fermions over 
|$o(A^))®|0)); we denote with |m)(g)|0)) these states. At weak 
coupling, we may truncate the Hilbert space to the simplest 
excitations, i.e., to states involving 1 e-h pair or 1 phonon cre- 
ated over the |7n)(8)|0)) states. We define the \m)(gi\q)) states, 
obtained by creating a phonon denoted by g = (r;, u) over the 
I'm) (g) |0)) states. Finally we introduce the \a) (g) |0)) states, 
obtained from the |r?i) (g) |0)) states by creating 1 electron- hole 
(e-h) pair. 

The approximation can be systematically improved by in- 
cluding two or more electron-hole and excitations in the trun- 
cated Hilbert space, at the cost of heavier computation. 

We now expand the interacting ground state in the trun- 
cated Hilbert space 

|*o(iV-F2))=^a™|m)®!0))-f 

m 

+^a™,,|m)(g)|g))-f ^ac!a)g)|0)) (7) 



and set up the Schrodinger equation 

H,i^i^u\^o{N + 2)) = £ |*o(iV + 2)). 



(8) 



We now consider the effects of the operators Hq and Vtot on 
thetermsof |*o(A^ + 2)). Choosing the |m)(g|0,q)), |a)(g|0)) 
states to be eigenstates of the noninteracting term Ho we have 



Ho|m)g)|0)) = £',„|m)(g>|0)), 
Ho\m)(g>\q)) = {Em+ujq)\m)(g)\q)), 
HQ\a}(g)\0)) = Ec\a)(g)\0)). 



(9) 

(10) 

(11) 



Let us consider the action of V and W on the same states, 
taking in account that V creates or annihilates up to 1 phonon 
and 1 e-h pair, and W is diagonal in the phonon states and 
can create or destroy up to 2 e-h pairs. 

{V + W)\m) ® |0)) = Y, V1„M) ® \l)) + 

m' ,q 

+ Y,^rn,n.'\m') g) |0)) + ^ W„,,^\a) (g) |0)) (12) 



(y + iy)|m)®|g)) = ^y^,^,|m')®10)) + 
+ Y,V:h,c.\a)(g>\0)) + Y,^n.,,n'\m')(g>\q)) (13) 

a m' 

{V + W)\a) ® |0)) = Y, Vl^,\m') ® \q)) + 

m' ,q 

+ Y W"<.,.'la') eg) |0)) + J2 Wo.,,^\m) g) |0)) . (14) 

a' m. 

The Schrodinger equation yields three coupled equations for 
the coefficients a's: 

{Err, - E)arrr + ^ »,„/ W„,,„/ + ^ Om'^qV;^!^^/ + 
m' 111'' ,<i 

-f^a„W^™,<,=0; (15) 

a 

{Em +U!q- E)a„r.q + ^J «m' .qW^™,™/ + 2_/ """' ^"".m' + 
m' m,' 

+ ^a.y,^,„=0; (16) 



{Ea - E)aa + ^ aa,Wa,a' + 
a.' 



(17) 



We define renormalized eigenenergies E'^ by taking a linear 
combination of the a states in such a way that 



{Ho + W)a,a' = Sac'Ea, 

Eq.(17) can be solved for the aa coefficients: 

pi p C/, amWa,m + 2_^ am,qV^^rrr) ■ 



Ea — E 



(18) 



(19) 



Sobstituing aa in Eq.(16) one gets: 

{Em +LOq- E)am,q + 2J "m' ,9 W^m.m' + 2j "™' ''^m.m' + 
m' m' 

v v' 

m' ,q' ol 

m' a 

Here, two important simplifications allow to proceed. First, 
as in Ref. [10], Eq.(49), one can show that 



Wm,m'=W^^^^,^5r^,r^,WF, 



(21) 



W I is the direct interaction among the particles forming 
the pair and Wf comes from the average over the occupied 
states on the Fermi sphere and is a 7n-independent constant; 
since W \^, vanishes for the W = Q property, it holds 



E' 



'm' ,q '* m.m 



,Wf. 



5 

m' ,q' 


'"'^2^ E'^~E '' 

a. 




\ ^ 1 •' m.a. \ 



(22) and therefore Eq.(26) must be solved self-consistently. Let us 
examine in detail the structure of the S[E] contribution. The 
matrix elements S^^rn' may be written as 



\y'^cii}rn,,m' I ^ m'^m.m' 



(29) 



(23) 



Indeed, V" is a one-body operator for the fermions; so the 
electron-hole pair in the a state must be created by one V 
factor and annihilated by the other; in this way, the W = Q 
pair is not touched. With these simplifications, the contri- 
butions in Eqs.(22) and (23) can be taken over to the l.h.s. 
of Eq.(20), where they just renormalize the eigen- energies of 
the \m) ® \q)) states. Thus, Em -\- uiq ^^ E'^, + (^q, and the 
Eq.(20) can easily solved for the am.q's, as we did for the a^'s 
in Eq.(19): 



where W^s is the true effective interaction between the elec- 
trons in the m states, while the other term represents the 
forward scattering amplitude F. 

The first-order self-energy WF5m,m' and the forward scat- 
tering term FmSm,m' are diagonal in the indices m and m' , 
and therefore they renormalize the non-interacting energy Em 
of the m states: 



E„ 



r(«) 



Em + WF+ Frr 



(30) 



E'„ 



E ^^ 






„' + 



E- 



EL-E 



If the effective interaction Weft is attractive and produces 
bound, localized states the spectrum of the Schrodinger equa- 
tion with the Hamiltonian in Eq.(27) contains discrete states 
below the unpaired states. In an extended system, we have 
bound states below the threshold of the continuum. The 
threshold may be defined (in clusters and in extended sys- 
tems) by 



(24) 



E, 



(R) 



T =mm[E\:^\E)], 

{m} 



(31) 



Finally, substituting Eqs.(19, 24) into Eq.(15), we can write 
the Schrodinger equation in terms of only the \m) states, with 
the excitations-mediated interactions and with renormalized 
quantities: 



= {Em - E)am + amWp + 



which, according to Eq.(30), takes into account all the pair- 
wise interactions except those between the particles in the 
pair. Note that this is an extensive quantity, i.e. an A'^ + 2- 
particle energy. The ground state energy E may be conve- 
niently written as 



+ 



E 



+2 > a„ 

m' ^vn" ,q,Oi 



E "" 

vn' ,vn" ,q 
m' ,ol 



VI 



,,v\, 



9 
in'' .vn' 



E = E. 



(R) 



-A; 



(32) 



E' „ +UJq- 
Wm,<yWcc.m' 



EL 



„V\, 



E 
W^ 



E 

+ 



{E'+Uq-E){E'^-E) 



A < indicates a Cooper-like instability of the normal Fermi 
liquid and its magnitude represents the binding energy of the 
pair. 

Below, we solve Eq.(26) explicitly for the CUO4 cluster with 
open boundary conditions, where the above theory is readily 
applied. 



E 



m' ,m" ,q,a,a' 



Wm.o.Vl^„VZ„^^,W^',m' 

{E'+LOq-E){EL-E){E'~E) 



(25) 



Here, the last two terms are of higher order and must be 
dropped; E is the ground state of the system with A'' -I- 2 
Fermions; yet, Eq.(25) is of the form of a Schrodinger equation 
with eigenvalue E for the added pair. We interpret am as 
the expansion coefficients over the W = Q pairs of the wave 
function of the dressed pair, j<^) = ^ am\'m) ® |0)). This 
obeys the Cooper-like equation 



-ffpairjlyS) = E\l^) 



(26) 



with the same E as in Eq.(8), but with an effective two-body 
Hamiltonian: 



H^ 



Ho + Wf + S[E] , 



(27) 



Here S is the U— dependent effective scattering operator 



5[S]. 



E'"' 



.W^ 



m,,af»'a,m' \ ^ 






m.m" m" ,m' 



E'a — E 



Z^ E' 



S' 



(28) 



IV. PROTOTYPE CLUSTER 

As an illustrative application of the above pairing scheme, 
in this preliminary work we focus on CUO4, the smallest clus- 
ter yielding W = Q pairing in the Hubbard model. This re- 
quires 4 holes, (total number, not referred to half filling); such 
a doping is somewhat unrealistic, but larger C4i,-symmetric 
clusters and the full Cu02 plane also show VK = pairing in 
the doping regime relevant for cuprates [9] [26] . Remarkably, 
in the pure Hubbard model, one can verify that A = A(4) at 
least at weak coupling [9] , which demonstrates that A has the 
physical meaning of an effective interaction. CUO4 represents 
a good test of the interplay between electronic and phononic 
pairing mechanisms since we can compare exact diagonaliza- 
tion results with the analytic approximations of the canoni- 
cal transformation. A further merit of this model is that it 
demonstrates dramatically the decisive role of symmetry in 
the electronic pairing mechanism: any serious distortion of 
the square symmerty restores the normal A > situation [9] . 
Since vibrations cause distortions it is not evident a priori 



that they tend to help pairing; in particular we may expect 
that Jahn- Teller distortions are going to prevent W = pair- 
ing altogether. On the other hand, the Frohlich mechanism of 
conventional superconductivity is based on phonon exchange. 
This suggests that the role of EP coupling is complex. 

Cu04 allows only the coupling to phonons at the centre or 
at the edge of the Brillouin Zone; however, phonons near the 
edge are precisely those most involved [14] [15]. Even in this 
small system the virtually exact diagonalizations are already 
hard and the next C4„-symmetric example, the CU5O4 cluster 
[27], is much more demanding for the number of vibrations 
and the size of the electronic Hilbert space. 
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FIG. 1. Pictorial representation of the ionic displacements in 
the eight normal modes of the Cu04 cluster, labelled according to 
the irreps of the C41, Group. 



Starting with the C4i,-symmetric arrangement, any dis- 
placement of the Oxygens in the plane can be analised in 
irreps, Ai, A2, Bi, B2, Ei, E2, see Fig.l. 

We suppose that the hopping integrals depend only on the 
bond lengths*. Hence, the EP coupling is expressed through 
just two parameters g and gox, defined as follows: denot- 
ing e.g. by t^ the hopping integral between Oxygen 1 and 
the Cu and by tl'^ the one between Oxygens 1 and 2 (see 



*Some Authors use an alternating sign convention for the 
bonds from a given Cu site. However, this is just a gauge; in 
the present CUO4 case, this corresponds to changing the sign 
of two opposite Oxygen orbitals. Even in the full plane, start- 
ing from positive t integrals, one can introduce staggered signs 
by negating a sublattice of O orbitals; then, one can arrange 
opposite signs for the bonds of each O by simply negating a 
sublattice of Cu. All this has no physical implications, and in 
our opinion does not help to visualize the real symmetry of 
the problem. 



Fig.l), g = [^]o and g,^ = [s|^]o. We take g < 
since a positive Cu-O hopping integral decreases as the Cu- 
O distance is increased. On the other hand gox > 0, since 
physically the 0-0 hopping integral has the opposite sign 
with respect the Cu-O one. Following Eqs.(3-6), the second- 
quantized electron-lattice Hamiltonian reads 



V 



^ fuj^blb,-! + f 2j('^CTPi<T + h.c) 

ia 

+ to. }_^{plp^ + la + h.c.) + UiY, "It'^u' + ^t"'"!") + 
i a i 

+ ^e.(&l + ^)^^; (33) 

V 

where p]^ and Pia are the hole creation and annihilation op- 
erators onto the oxygen i = 1, .., 4 with spin a =1, 1, dl and 
da are the hole creation and annihilation operators onto the 
central copper site, while n'^ — pl^Pia and ria — d^da 
are the corresponding number operator. Henceforth we set 
Ep = Ed ~ for convenience, since this simple choice is ade- 
quate for the present qualitative purposes. Also, we are as- 
suming for simplicity that the Oxygen-Oxygen hopping tox 
is zero, and 0-0 hoppings are important only once the ions 
are moved. Similar results are obtained using a realistic tox, 
except that pair binding energies are somewhat reduced. 
The Hn matrices are given by 

Hai = -g '^{dtpia + h.c.) + 

-Jk9ox ^(p^Pi + la + h.c.) ■ 

i a 

Hbi = -g 2_^id'ipi<T - dlp2a + dlpscr ~ dlp4a + h.c.) ; 

a 

Hb.2 = -7^9ox ^(pLp2ct -pLp3<t + 

plaPi" -pLpi<^ + h.c.) ; 
Hei^ = —;=g y^^{-dipiadip3a + h.c.) + 

-gox 2J (-pLp2<t + pLp3<t + pLp4<t - pLpi<^ + h.c.) ; 

a 
HE^y = -J=9 ^^{dlp2a - dipAa + k.C.) + 



r,9ox 






^/2 



(pLp2<t + pLp3<t - pLp4<t - pIoPit + h.c' 



+pLpit + h.c.) ■ 

HE2y = ^gox ^(-pLp2<t -pLP3a 

a 

+plaP4a + p\oPl<y + h.c.) . 



(34) 



In order to make contact with the physics of cuprates, let us 
discuss the connection between the normal modes of the Cu04 
cluster andthe phonon modes of the Cu-O planes. There is 
experimental evidence [15] that the possibly relevant modes 
for superconductivity lie on the Cu02 planes and have a Cu- 
O bond stretching origin. In particular the LO half-breathing 
mode with k — (vr, 0), (0, tv) is believed to couple significantly 
with the doped holes in the superconducting regime. In the 
Cu04 cluster the half breathing modes are contained in the 
breathing mode Ai and in the quadrupolar mode Bi by means 
of the linear combination qa^ ± gs., . We argue that qualita- 
tively the effect of the coupling with the Ai and Bi modes 
should give us clues about the interplay between electronic 
W = pairing and phonon exchange. 



V. LOWEST-ORDER EFFECTIVE 
INTERACTION IN CUO4 

The mere Hubbard CUO4 cluster with 0-0 hopping tox ~ 
yields [9] A(4) < 0, due to a couple of degenerate W — 
bound pairs, in the Ai and B2 irreps of the C^v group; 
therefore in Eq.(28) we set the m = m' labels accordingly. 
At weak coupling, we may simplify Eq.(28), neglecting all 
renormalizations; the phonon-mediated interaction for the B2 
pair reads: 



^k 



m," ,q 



,K 



m" ,B2 



+ iu'q 



E 



= -4go^ 



2eAi +^S2 



E 



(35) 



Note that in the denominator in the r.h.s., 2eAi + ^B2 is the 
energy of an unrenormalized excited |m)®|g)) state, which at 
weak couplig is higher than the ground state energy E; hence 
the r.h.s. must be negative and the B2 phonon is synergic with 
electronic pairing. On the other hand, the vibronic effective 
interaction for the Ai pair is: 



E 



2A1 



2eAi + ^.4i 



yq yq 

Ai,m" m" ,Ai 


4 2 


( ^k 




Em" +oj'q-E 


\2£Ai +(^B2 


-E 


- E 2eAi + ^Ei 


-E 


2eAi +^E2 


- E 



+ 



This shows that in the Ai sector the total sign depends on 
the relative weight of attractive and repulsive contributions. 
Eqs. (35,36) show that at weak coupling Ai and B2 modes are 
synergic to the W = pairing, while both longitudinal and 
transverse E modes are pair-breaking. The half-breathing 
modes that are deemed most important [15] [16] are Aj ± Bj 
combinations, but Bi does not appear in Eqs. (35, 36). The 
numerical calculations reported below confirm these findings 
in a broad range of parameters. 

For the sake of argument, in the explicit calculations we 
took all the normal modes with the same energy eo = huo = 
10~^eV and A,, = 1. This sets the length scale of lattice effects 

"^ A where we used M = 2.7 x 10"' 



Co = 



10" 



2i\/u 

for Oxygen. 

With this choice, the Cooper-like equation (25) reads 



^Kg 



(2eAi ~E)- 



16 



esi + £Ai 
1 



E 



4 2 

3^°" 2eAi +^o-E 



in the Ai channel and 



(2£Ai - E) 



16 



4,gL 



1 1 

2 £Ai + £a[ 



1 1 



E 



(37) 



£Bi + £Ai 
1 



2£Ai +LOa - E 



E 2eAi 
= 0. 



- £a' - E 



(38) 



for B2 pairs. The eigenvalue E, like in Eq.(8), is the total 
energy of the cluster; it must be compared with the threshold 



r{«) 



(R) 



E)p of Eq.(31), whose noniteracting limit is _E-j, = 2eAi 
since the degenerate level energy (see Table II, Appendix A) 
is £p — 0. It turns out that using Eq.(32) in the weak coupling 
approximation, that ignores renormalizations, the effective in- 
teraction is A = _E — 2eAi i i^i Appendix B we verify by pertur- 
bation theory that like in the Hubbard model, A = A(4); this 
supports our interpretation of A as minus the pairing energy. 
The trend of A in both channels is shown in Fig. 2. The 
vibrations split the degeneracy of the W = pairs, effectively 
lowering the symmetry like a nonvanishing tox- Pairing is 
enhanced in the Ai sector as well, albeit less than in B2; 
without phonons, A ~ — 20meV for both the W = pairs. 
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FIG. 2. Analytical results of the canonical transformation: 
pair binding energy in the Ai and B2 sectors as a function of Qox- 



Here we used A,, = 1 for every mode, t = leV, tox 

0-1 
gox is in units of eo/Co = IcVx A , A is in eV. 



■ 0, U = leV; 



VI. JAHN-TELLER MIXING OF ELECTRONIC 
GROUND STATES AND PAIRING 

If really the inclusion of the lattice degrees of freedom sys- 
tematically leads to A > 0, purely electronic cluster models 
become totally irrelevant to superconductivity, as it was ar- 
gued [18]. However small, the CUO4 cluster yields electronic 
pairing and allows to test this important point. 

In this Section we set up a conventional calculation of the 
JT effect involving degenerate electronic ground states and 
their mixing with the vibrations. We first take the nuclei as 



frozen in a C4„-syninietric configuration and diagonalizc the 
purely electronic part of the Hamiltonian: 



^e?°'=*^(4p-+ft.c) + 



^(E 



-if-if+^r-D 



(39) 



7-CUO4 



As before, we are using tax = Ep = £d = 0- 

The JT effect arises if the ground state of ^^"^'* is degener- 
ate, with a ground state multiplet {|*I'i), • . • , |*I'n)} such that 
Hfi^'^-'l-^k) =Eo\^k}- If we take matrix elements of H^j^^^u 
in this truncated (n- dimensional) electronic basis integrating 
over electrons and keeping boson operators we get the dynam- 
ical JT Hamiltonian, [29], [30] with matrix elements: 



Ha'!p = {Eo + ^ hojy^blb,-,) Sa.0 + (*«| V'l'f/s) , 



(40) 



It is worth noting that neglecting the nuclear kinetic energy 



(i 



2 A/ Z-ti 



) and treating the nuclear positions as 



variational parameters corresponds to the static JT Hamilto- 
nian, but we follow the dynamic treatment which is superior. 
In the following we assume that the initial configuration is 
stable with respect to the mode Ai which only changes the 
scale of the CUO4 molecule. Since this mode does not pro- 
duce any JT distorsion, it is not involved in the arguments of 
Ref. [18]. In this Section, we study A(4) in this approxima- 
tion, according to Eq.(l). The ground state with 2 holes is 
a nondegenerate totalsymmetric singlet unaffected by the JT 
effect; in the other cases, the use of the Hamiltonian (40) is 
justified provided that the excited states are several phonon 
energies above the ground state. 



A. Three-hole ground state mixing 

With three holes the ground state belongs to the 3- 
dimensional irrep of S4 which in C4t, breaks into Bi (B E. 
To illustrate the electronic structure and its dependence on 
distortions, in Fig. 3 we show the adiabatic potential energy 
surface projected along the B2 distortion. Projecting on the 
other directions, we obtain similar trends. It is clear that 
the ground state multiplet is well separated from the excited 
states and hence this treatment of the JT effect is well justi- 
fied. 




-2.95 \ 
FIG. 3. Adiabatic potential energy surfaces along B2 for the 
ground and first excited 3-hole states of CuOa- Here, gs^ denotes 
the classical normal coordinate; U = 5eV, t = leV, g = —2.4, 

o 

gox = 0.6 in units eV/A and a;,, = O.leV V/y; the dispaccment qsj 

o 

is in A, energies are in eV. The ground state multiplet is below tiie 
first excited state by ~ leV. 



Since E ® E contains all the irreps of dv, all the normal 
modes are JT active in this case. Following Eqs.(40), we com- 
puted the following V matrix elements in the 3-hole ground 
state multiplet with the Hubbard interaction, using Eq.(34). 
The 4 independent elements at the optimal value U /t ~ 5, 
where the W = G pair binding energy is maximum, are: 



71 = (*Sll^i5lJ*i5j = 0.17 <7o, 

72 = {*Si1^£;2J*bJ= 0.245 + 0.17,90. 

73 = (*Bj-ffi3il*i5j =0.243 

74 = (^-Bj-ffsal^isJ = -l-Q^Qox 
The JT Hamiltonian reads 



(41) 
(42) 
(43) 
(44) 



Hjt(3) = 



Eoii) + ^nu^b% 



l3x3 + 



+ ^in{bl + hr,)M^ 



(45) 



where Eo{'i) is the ground state energy of H^^ ^ with 3 holes 
and 



/ 
Mbi = 73 1 
V -1 



(46) 



/ 

Mb2 = 74 1 

V 1 



(47) 



Me,^ 



Me,, 




(48) 



(49) 



and 



/ -11 
Me^, =71-1 

V 1 



/oil 

Me^, = 71 1 

V 1 



(50) 



(51) 



We numerically diagonalized Hjt{3) in the Hilbert space 
spanned by I'^-q) ® |$jv h), where rj = Bi, Ex, Ey is the elec- 
tronic state and |$jv ^) is a vibration state in the truncated 
Hilbert space with all modes having vibrational quantum 
numbers < Nph- Excluding the breathing mode the size of 
the problem is S{Nph + 1)^. We consider Nph = 3, since 
already in the weak coupling regime A changes sign. We 
studied a;,, = a; = O.leV^ A,, = 1 V77 in the range of \g\ and 

o 

\gox\ between and 2.4eV/A, which means that the number 
ratios llilyj^j^/fi-'^ vary between and 0.5. This condition 
ensures that the weak coupling regime behooves and Nph ~ 3 
is indeed adequate. The results are shown in Fig. 4. 




FIG. 4. (a) Vibronic correction AE to Eo{3) to the E state 

(lower surface) and Bi state (upper one) as a function of g and 

gox- (b) the same surfaces from another point of view. Here, 

t = leV, U = 5eV, AE is eV of hw, g and gox are in units of 

0-1 
£o/?o = leVx A . 



The lower surface represents the ground state energy shift 
AE for electronic states belonging to the degenerate irrep 
E, while the higher one is AE refers to the Bi state. The 
JT effect partially removes the three-fold degeneracy and the 
ground state is a -E doublet. Note that according to the text- 
book, static JT effect, one should observe a total removal of 
the degeneracy. This is however not borne out by the dynam- 
ical calculation and, for gox = 0, all the three states remain 
degenerate (see Fig.4.b). 

The way the system dynamically distorts is also of interest. 
The only vibration having the coordinate on the diagonal of 
Hjt{?>) is B\; thus, the E doublet can only distort along 
the Bi normal mode. In other terms, with qsi = Csi (&b + 

&flj, (gBi) = (*|jqBj*k) = -(^Ijgsjn,} ^ 0. The 
trend of (gsi) as a function of g and gox is shown in Fig. 5. 
We observe that ((Jsi) ^ as (/ — > 0; we also remark that 
the deformation depends essentially by g and only weakly on 
gox because according to Eq.(43) the coupUng constant 73 
responsible for the distorsion along Bi depends on g and not 
on gox- The E-B\ splitting, on the contrary, depends on gox 
and only weakly on g. The naive expectation that splittings 
go along with distortions only holds for static ones. 




1.44 



FIG. 5. {qBi) as a function of —g and gox- Here t = leV, 
U = 5eV, {(JBi) is in units of ^o, g and gox are in units of 

0-1 
£o/Co = leVx A 



B. Four- hole ground state mixing 

With 4 holes, -ff^j" *" has a two-fold degenerate ground 
state; it belongs to the two dimensional irrep of S4 that breaks 
in Ai(BB2 in C41,. Thus, the only JT-active mode is B2, which 
makes the problem exactly resoluble in terms of a continued 
fraction [31]. 

Following again Eqs.(40) the second-quantized JT Hamil- 
tonian with four particles reads 

Hjt{4.) = [hujB2 bl^be-, + So (4)] (g) l2x2 + 

1 



Cfl2 75(&l3, +bB2) 



1 



(52) 



in the space spanned by the electronic ground states. The 
coupling constant 75 is given by 



-75 = {'i' A^\H B^l^B:,} = 1-19 gox 



(53) 



where as usual, the matrix element in Eq.(53) is evaluated at 
U/t ~ 5. 

The ground state energy of Hjt{4:) in the sector of sym- 
metry 77 coincides with the lowest pole of the Green function 



G...(S) = {{0|®(*,| 



1 



E-Hjt{'1)+iO+ 
1 



l*.)®|0)) 



x(o) 



(€b2^5)2 



(54) 



x(i)- 



HIb-.-'b)^ 



5(iB2^!>r 



x(4)-... 



where rj = Ai,B2 and x{i^) = E — -Bo(4) — nhuiBz- G^,^(i5) 
does not depend by ri and the energy corrections in the Ai 
and B2 sectors are the same. The ground state energy shift 
AE as a function of gox is plotted in Fig. 6. 



AE 




FIG. 6. Ground state energy shift AE as a function of gox 
Here t = leV, U = 5eV, A_B is in eV, gox is in units of 



eo/?o = leVx A , A^ 



1, f^s. 



: O.leV. 



FIG. 8. A as a function of —g and gox, according to the the- 
ory of the present Section. A is in eV, g and gox are in units of 

o-l 

£o/Co = leVx A Here t = leV, U = 5eV, ojb^ = O.leV. 



There are no diagonal couplings in Eq.(52), which implies 
no distorsions; the energy correction is much smaller than in 
the three holes case. 

In Fig. 7 we show the adiabatic potential curve, with the 
same parameters as in Fig. 7. In contrast with the three-hole 
case, the ground state multiplet is separated from the excited 
states by ~ 100 -r- 150meV which is comparable the phonon 
energies. Hence we expect that in this case the approxima- 
tion restricting the Hilbert space to the lowest multiplet is 
not justified. However, in the case of weak coupling to soft 
vibrations, with uj small compared to the gap in the electronic 
spectrum, this approximation should work well. 
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FIG. 7. Adiabatic potential energy curve for the low lying 
4-hole states of CuOa. Here, qg^ denotes the classical normal 

o 

coordinate; U = 5eV, t = leV, g = —2.4, gox = 0.6 eV/A and 

° "^ ° 

^Mujg = 20eV/A ; the displacement is in A, energies are in 

eV. The ground state multiplet is below the first excited state by 

~ 100 -=- ISOmcV which is of the order of 2|A|. 



Both the systems with four and three holes gain energy 
by the JT effect; A remains negative only in the weak 
EP coupling regime, since the decrease of 2E{3) overcomes 
the decrease of -E(4). This is due to the fact that the 
system with three holes can gain energy by mixing with 
Bi,B2,Eix,Eiy,E2x,E2y vibrations, while the system with 
four holes can do it only with the B2 mode. Moreover the fac- 
tor 2 in front of E{3) in the expression for A further favours 
the distorsion in the 3-holes case. 

These results in line with Ref. [18], would imply that the 
electronic pairing is limited to relatively weak EP couplings 
and that at any rate the vibrations are invariably detrimental 
to W = pairing. However we know from Sect.V that this 
conclusion is remarkably but definitely wrong, because the 
full theory predicts synergy of vibrations and W = pairing 
at least at weak coupling. This failure of the JT Hamiltonian 
is due to the neglect of the electronic excited states, causing a 
severe overestimate of the 4 hole energy. The physical reason 
is that if we restrict the mixing to the degenerate states the 
electronic wave function is too rigid. On the other hand, in- 
cluding all the multiplet of states arising from the degenerate 
one-electron level the pair can achieve the fiexibility which 
allows it to follow adiabatically the vibration-induced defor- 
mations, as we shall see in the next Section. 

Detecting the failure of a textbook procedure is by itself a 
potentially very interesting result; it is a merit of a relatively 
simple model like this that allows to understand in detail how 
this arises. 



VII. NUMERICAL RESULTS OF THE FULL 
THEORY 



C W = pairing in the presence of Jahn- Teller 
distortions 

Collecting together the results of the present Section we 
obtain the behaviour of A in a popular approximation that 
neglects the excited electronic states. In Fig.8 we show the 
plot of A as a function of —g and Qox, shown with the zero- 
energy plane. The maximum distorsion along Bi compatible 

with A < (see Fig.S.b) is (g™"^) ~ 3 x 10"^ A, which is 



attained at 17 ~ leVx A and gox 



0. 




Since the conventional JT Hamiltonian is not enough, to 
see what really happens in this model with increasing EP 
coupling, where the analytic treatment of Sect.V loses va- 
lidity, we resort to numerical methods. In this Section we 
explore the pairing scenario numerically, which offers an inde- 
pendent check of the weak-coupling calculations and permits 
to go beyond the weak coupling regime. First, we analyze 
one phonon at a time (A,, — 1), turning off the all others 
(A,, — 0); Lu^ = 0.1 for all modes. In this way we see which 
kind of phonon is cooperative with the W = pairing and 
which is not. We have performed these calculations in a vir- 
tually exact way, by including a number of phonons Nph up 
to 20. To this end we take advantage of the recently proposed 
spin-disentangled diagonalization technique [11]. The results 
are shown in Figs. 9, 10 and 11 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 

FIG. 9. Exact diagonalization results for A(4) in cV, with only 
the Ai phonon active, Npf^ = 20, as a function of gox for different 
values of g: g = —0.2 (crosses); g = —0.5 (triangles); g = —1 
(diamonds). Here we used t = leV, tox = 0, U = leV, gox and g 

o-l 

are in units of eo/Co = leVx A 
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FIG. 11. Exact diagonalization results for A(4) in eV, with only 
the El phonon active, Np/^ = 20, as a function of gox for different 
values of g: g = —0.2 (crosses); g = —0.5 (triangles); g = —1 
(diamonds). Here we used t = leV, tax = 0, U = leV, gox and g 

o-l 

are in units of £0/^0 = leVx A 




oox 

FIG. 10. Exact diagonalization results for A(4) in eV, with 



N, 



ph 



20 



only one phonon active, ^2 (triangles) or E2 (diamonds), 

. Here we used t = leV, U = leV; gox and g are in units of 

0-1 
£o/Co = leVx A 



The plots show the trend of A(4) as a function of gox and g. 
It appears (see Figs. 9 and 10) that if gox is increased, the Ai 
and B2 phonons enhances the pairing , even beyond the weak 
coupling regime. The further enhancement of |A(4)| due to 
Ai as \g\ is increased is not predicted by the weak coupling 
theory [Eq.(36)]. The Bi phonon is slightly suppressive, but 
it affects the pairing energy on a scale of 10~^eV and hence 
its contribution is negligible. On the other hand the (longi- 
tudinal and transverse) E phonons (see Figs. 10 and 11) have 
an unambiguos tendency to destroy the pairing. In particular 
the El mode does it both by increasing gox and by increasing 

The presented results show that the behaviour of the in- 
dividual phonons is essentially the same as predicted ana- 
lytically in Sect.V: some of them (^1 and B2) act in a co- 
operative way with the electronic pairing mechanism; some 
other {E phonons) does not; the Bi mode is quite inopera- 
tive. However, for a proper understanding of the conflicting 
vibronic effects we need to include as many phonon modes as 
possible at the same time. In this case the exact diagonaliza- 
tions become hard even with a modest Nph per vibration. We 



performed exact diagonalizations of H^ 



-latt 



with five active 



modes; with 4 holes, the size of the problem is lQQ{Nph + 1)^; 
we could afford Nph = 3 for each. Some results are shown in 
Figs. 12. 

In Fig. 12. a we included the vibrations with rj = 
Ai,Bi,B2,E2x,E2y; one notes a strong, monotonic increase 
of the binding energy with both gox and \g\. The weak- 
coupling theory of Sect.V qualitatively explains the gox de- 
pendence but not the \g\ one: when the EP coupling gets 
strong, the Cu-O stretching grows important. In Fig.lO (dia- 
monds) we noted that the E2 vibrations alone tend to destroy 
pairing; here we observe that when they compete with Ai and 
B2 their effects are utterly suppressed. It is possible that the 
couplings to the pair-breaking E modes are somewhat under- 
estimated by the choiche of A parameters. 

In Fig. 12. b we included the vibrations with 77 = 
Ai,Bi,B2,Eix,Eiy, and we observe that A now becomes 



10 



positive at moderate Qox- Comparing with the above results 

on individual modes, we observe that in going from Fig. 12. a 

to Fig. 12. b we are replacing the pair-breaking transversal 

E2 phonons by the pair-breaking, longitudinal Ei modes. 

We conclude that the longitudinal ones are more efRcient 

in restoring the repulsion and at intermediate coupling they 

overwhelm the pair-healing A\ and 52- 

0-1 
However for Qox ^ 0.15 eV x A , the cooperative modes 

win and A (4) gets negative again. This remarkable behavior 

could not be anticipated by the weak coupling approach of 

Sect.V, where only the one-phonon exchange diagrams were 

included as in the BCS theory. 

In Fig.l2.c the active modes are 77 = Bi,Eix,Eiy, E2x, E2y\ 
these are all pair-breaking individually and switching them all 
together, they readily unbind pairs, leading to strong positive 
A. However, unaspectedly, we again find attraction at large 
enough Qox- 

It is likely that the couplings to the pair-breaking E modes 
are somewhat overestimated by the choiche of A parameters 
in Fig.l2.b,c. The pairing at strong coupling observed in 
Figs.l2.b,c results from more complicated interactions leading 
to bipolaron formation. This recalls the charge-ordered super- 
lattice phase found in Ref [25] ; however they used a Hubbard- 
Holstein model and, since the system is one-dimensional, the 
electronic pairing does not occour in their case. 
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FIG. 12. A(4) in cV as a function of gox for different values 

of g. A,, = 1 all the vibrations, exeept: A,, = for rj = Ei (a); 

A,, = for ri = E2 (b); A,, = for 77 = Ai, B2 (c). Npi^=3 for each 

active mode; g = —0.2 (diamonds); g = —0.5 (triangles); g = —1 

(crosses). Here we used t = leV, tox = 0, U = leV, gox and g are 

0-1 
in units of eo/?o = leVx A 
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VIII. CONCLUSIONS 

Introducing vibrations and vibronic couplings into a 
strongly correlated model opens up a rich scenario where, 
among other possibilities, pairing can be achieved by a syn- 
ergy of electronic correlation and phonon-exchange. A possi- 
ble outcome, however, is competition among different symme- 
try vibration modes and electronic excitations. We illustrate 
the situation by using a Cu04 model that allows a full treat- 
ment of all degrees of freedom and hosts bound W = pairs 
when undistorted, has vibrations of the same symmetries as 
the CuO plane and is numerically affordable. A popular recipe 
for computing JT distorted molecules prescribes restricting 
to the degenerate electronic levels letting them interact with 
the JT active modes. A static treatment invariably leads to a 
complete removal of the symmetry and a nodegenerate ground 
state. We put forth a fuller dynamical theory which partly 
preserves the degeneracy; however, the vibrations are always 
opposing W = pairing which is thereby reduced to a weak 
EP coupling effect. This restricted basis, however, may only 
be valid provided that the excited states of the unperturbed 
electronic Hamiltonian are far removed from the ground state 
on the energy scale set by the frequency of the relevant phonon 
modes. With the cuprates in mind, we consider a situation 
when the phonon energies and the superconducting gap are 
comparable, in the 0.1 eV range; we diagonalize the full model 
keeping up to 5 simultaneous modes and vibrational quantum 
numbers Nph < 3. Depending on the parameters, a rich phe- 
nomenology emerges from the numerical data. Pairing pre- 
vails at weak EP coupling, as expected, but the phonon contri- 
butions which dominate in such a case turns out to contribute 
to the pairing rather than opposing it. The correct trend is 
predicted by a canonical transformation approach, which also 
explains the pairing or pair-breaking character of the modes. 
In particular it is found that the half-breathing modes give 
a synergic contribution to the purely electronic pairing; since 
they are believed to be mainly involved in optimally doped 
cuprates, our findings suggest a joint mechanism, with the 
Hubbard model that captures a crucial part of the physics. 

This agreement validates the canonical transformation ap- 
proach, which allows to carry out useful calculations even in 
large systems that o not lend themselves to exact diagonal- 
ization. 

At intermediate coupling the outcome of the theory de- 
pends essentially on the relative weight of the coupling to the 
longitudinal and transverse vector modes, which destabilize 
pairing most effectively. 

Remarkably, however, the vibrations restore pairing again 
at strong coupling, when a bipolaronic regime prevails. This 
scenario was also drawn in the context of an extended t-J 
model, where the half-breathing mode was found to enhance 
electronic pairing [17] [32]. 

Finally, experimental data on nanopowders [33] also indi- 
cate that one should not be overly pessimistic about cluster 
calculations. The pairing that shows up there can be relevant 
and physically insightful concerning the interplay of various 
degrees of freedom on pair structure and formation. 



IX. APPENDIX A:W = PAIRS IN THE CUO4 
CLUSTER 

The Cu04 Hubbard Hamiltonian has C4„ symmetry. When 
the Oxygen-Oxygen hopping is absent, the symmetry group 
is the permutation group 5*4 , and although for convenience we 
continue to use the subgroup C4„ labels, it is S4 that must be 
used for the W = theorem [12]. The character table reads: 



C4i, 




C2 


r(+) r(-> 


O-x, CTy 


a+, (T_ 


Symmetry 


Ai 




1 


1 


1 


1 


x'+y' 


A2 




1 


1 


-1 


-1 


(x/y) - {y/x) 


Bi 




1 


-1 


1 


-1 


x'-y' 


B2 




1 


-1 


-1 


1 


xy 


E 


2 


-2 











{x,y) 



Table I: Character table of the C4v symmetry group. Here 1 



the 
the 



denotes the identity, C2 the 180 degrees rotation, C\ , C\ 
counterclockwise and clockwise 90 degrees rotations, a^, cry 
reflection with respect to the y = and x = axis and o-^-, u- the 
reflection with respect to the x = y and x = —y diagonals. In the 
last column we show typical basis Junctions. 



Setting for simplicity Ed = £p = 0, the one-body spectrum 
of the Cu04 Hamiltonian has the following eigenvalues: 





t 


t 


£E„ 
t 


t 


r - V4 + r2 


-2r 








r + V4 + r2 


-2 











2 



Table II: One-body levels of the CuO^ cluster in units of t; the 
last line reports the values for tox = which are used in the text. 



Here we label the eigenvalues by the irreps of the corre- 
sponding eigenfunctions. The level energies are in units of t, 
T = tox/t; for tox = 0, ee^ = £b„ = ebi = 0. The corre- 
sponding one body creation operators of a particle in each of 
these eigenstates are: 



-Eya ^ 



(pL 



(pL 



-via 



pI 



via 



-P\a 



4.(2) = 



-r^ {adi +pI +pL+pL+pL) 

(M+pL+pL+pL+pL) 



(3^+4 



where a and /3 depend on r as follows: 
4 (-l-r^ + rV4-fT2) 



-5 T - 2 r3 -h V4 + r2 + 2t^ V4 + r^ 



(55) 

(56) 

(57) 

(58) 

(59) 
(60) 

(61) 
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/3 = 



4 (l + T2 + rV4 + r2 



5 T + 2 r3 + V4 + t2 + 2 t2 V4 + r^ 



(62) 



By the W = theorem [12], the irrep Ai(BB2 of the group 5*4 
which is not represented in the one-body spectrum must yield 
singlet eigenstates with no double occupation. Projecting one 
finds: 



At Let c1" + — 



,t J 



/g y^E^t'^Eyl +<^Ey]<^E, 



These are readily verified to create W = pairs. 



(63) 



A(4) 



Hi 

16 



1 



2(-eAi + £a' 



2 ^ 



(67) 



This must be compared with Eq.(38), Sect.V, that can be 
solved iteratively for A(4) = E — 2eAi inserting ebi ~ 
from Appendix A. The second iteration yields Eq.(67), sup- 
porting the identification A(4) — A(4) at this order. Both 
quantities represent the effective interaction of the dressed 
B2 W = pair. Indeed, much information about the ground 
state energies cancels out if one applies Eq.(l); the canonical 
transformation is a much more practical way to represent the 
efi'ective interaction. 



X. APPENDIX B: PAIR BINDING ENERGY 

Here we calculate A of the B2 pair in the CUO4 cluster by 
second-order perturbation theory in both W and V and com- 
pare with A obtained by solving Eqs. (37,38). Basically the 
same holds for the Ai pair. A^ = 1 throughout this Appendix. 

The ground state with two particles belongs to Ai and, 
using the notations of Table III, its energy reads 
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g^{d2 + 2d-i) - pL (da + di) + 2V2go^gd3. (64) 


[2 
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i = -EAi + ea; 


^ = -EAi +u;b^ 


J^ = -£-4l + ^E, 






^ = -£Ai + UJE^ 


^ = £.4', + ^E^ 


^ = Ea', + ^B2 





Table III: Shorthand notations used in Eqs. (64,65,66). 
With 3 particles the ground state is a _E doublet, and 



E{S)^2eA^ + ^U + U'' 
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2 1 1 
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(JJB2 2iOEi 2lue2 
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(65) 



The ground state with 4 particles belongs to B2, as predicted 
by the canonical tranformation; 
one gets: 



S(4)=2eAi + ^C/ + f/' 
16 



25 



29 



di 



128 EAi 512 

g'^ [d2 + d3 + ds] - 
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8 1 1 d3 d4 ds dfi 

^ \ h— ■ + — 4- — -h — 

UJB2 ^E, (j->Er, 2 2 2 2 



4- 



goxgV2 [ds - d2] . 
Finally using Eq.(l) and setting cj^ = tJo, we obtain 



(66) 
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